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Theoretical Descriptions of the Nucleon Electromagnetic Form Factors in 
the Space- and Time-like regions 
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Some recent advances in the theoretical description of the Nucleon electromagnetic form factors, both in the 
space- and time-like regions, will be briefly illustrated. In particular, both the present stage of the lattice calcula- 
tions and updated phenomenological approaches, like the ones based on dispersion relations and on microscopical 
models, will be reported. 



1. Introduction 

The Nucleon electromagnetic form factors (em 
ff's) represent a relevant chapter in our path to- 
wards the understanding of hadron dynamics. It 
is striking that, after almost half a century of 
experimental investigations, this field could still 
hold a number of surprises (e.g., the intriguing 
fall-off of the proton charge ff more rapid than 
the assumed dipole one), or have some kincmat- 
ical regions not fully explored (e.g., the timelike 
region, in order to obtain separate information on 
ff's ). In particular, let us remind (see, e.g., |H2j 
for recent reviews) that in the spacelike (SL) re- 
gion the existence of a possible zero in G V E for 
-q 2 = Q 2 > 7 GeV 2 (q 2 = the square mo- 
mentum transfer), as suggested by a linear ex- 
trapolation of the data obtained by the polariza- 
tion transfer experiments, could open new per- 
spectives. Indeed, the dipole form of the Nu- 
cleon ff's is supported by the pQCD, therefore 
a sharp difference between the charge distribu- 
tion of the current quarks and the charge dis- 
tribution of the "constituent" quarks appears a 
clear signature of NpQCD effects in act. Even 
low-Q 2 ff's (cf the forthcoming analysis of E08- 
007 at TJLAB) could give new stringent con- 
straints for builders of phenomenological mod- 
els. In the the timelike (TL) region, analogously 
relevant issues can be met. For instance, the 
neutron data around q 2 ~ 4 GeV 2 show a be- 



havior, \G 



eff\ 
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sharply different from 



the one suggested by pQCD, naively indicating 
\Geff\ ~ \ e d/e u \ \G p eff \. Moreover, TL proton 
data show many interesting structures (see, e.g. 
[3]), that could give a lot of information on the 
hadronic component of the photon wave function. 
On the theoretical side, in a very schematic 
way, one can list the following approaches: i) Lat- 
tice calculations, based on a Euclidean 4D space, 
that can address the SL region only; ii) analytic 
approaches, based on dispersion relations, that 
can investigate both SL and TL regions, even 
yielding predictions for the unphysical interval 
< q 2 < 4Mjy-; hi) purely phenomenological 
methods, with constituent quarks (CQ) as ba- 
sic ingredients, that can be exploited for describ- 
ing both SL and TL regions; iv) field theoretical 
approaches, like the ones based on the Dyson- 
Schwinger equation, that can be applied to SL 
and TL regions. In this short review we will 
restrict the presentation to i), ii) and iii) (see, 
e.g.,[T] for a different discussion from i) to iv)). 
It should be also mentioned a quite recent ap- 
proach based on the gravity dual models of QCD 
(or holographic QCD), like the one in Ref. [3]. 

2. Lattice "Data" 

In principle, Lattice QCD allows one to achieve 
a non perturbative description of the Nucleon ff's, 
as shown by the more and more accurate results 
obtained in recent years. In particular, such inter- 
esting outcomes have taken profit from: i) an im- 



pressive growth of the available computing power, 
that has allowed in nowadays results, e.g., to re- 
lease the quenched approximation (where the sea- 
quark loop contributions are disregarded); ii) new 
algorithms, leading to a better description of the 
fundamental symmetries of QCD, and last, but 
not least, hi) substantial improvements of the chi- 
ral extrapolation techniques. 

On the other hand, we should mention the main 
obstacles to be overcome: 

• in a 4D Euclidean Lattice (i.e., it — > t), 
one cannot directly evaluate the matrix el- 
ements of the em current operator. There- 
fore a careful study of the proper ratio be- 
tween three- and two-point functions is nec- 
essary in order to get " large values" of the 
Euclidean time, that eventually allows to 
filter the contribution from hadronic states 
with Nucleon quantum numbers. This fil- 
tering prevents the study of TL ff's, since an 
exponential fall-off in time lacks for q 2 > 0; 

• present Lattice data are obtained with val- 
ues of the coupling constant j3 = 6/g 2 and 
the hopping parameter, k, that lead to siz- 
able u, d quark masses. Usually, the calcula- 
tions are labeled by the corresponding pion 
mass, from which one can infer the quark 
mass in the actual data. In order to move 
toward the current quark masses, one has 
to reduce the Lattice spacing (that allows 
to accommodate lighter quarks) or to ex- 
ploit extrapolation formulas, obtained from 
the low-energy regime of the QCD. Unfor- 
tunately such chiral extrapolations have a 
well-known "non analytic" (log m^) behav- 
ior for small value of the pion mass. The 
lowest value of the pion mass in the present 
Lattice calculations for the Nucleon ff's is 
- 360 MeV [5]; 

• the smallness of momentum transfers in- 
vestigated in a Lattice is related to the 
space extent, L, namely \q\ ~ 2tt/L. 
This has influence on the range of appli- 
cability of the chiral extrapolations, about 
Q 2 < 0.4 GeV 2 , not always reached by the 
present Lattices (e.g., Q 2 ow ~ 0.23 GeV 2 for 
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Figure 1. LHP Coll. results for the "isovec- 
tor" (see text) approximation of the Dirac ff 
F\{Q 2 ) (upper panel) and for the "isovector" ra- 
dius (lower panel). Solid and dashed lines: differ- 
ent extrapolation formulas. (After Ref. [5])- 
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the LHP Coll. [5] and 

for the QCDSF Coll. [Sj and 

0.15 GeV 2 , for the Cyprus-MIT Coir [7], 

but with quenching approximation); 

• "large" momentum transfers in the range 
of few GeV 2 are in principle accessible 
by nowadays Lattice spacing, but Fourier 
transforms of two- and three-point func- 
tions become noise-dominated for square 
momentum transfers beyond 2 GeV 2 ; 

As an example of an updated Lattice calcu- 
lation of the Nucleon ff's, the results obtained 



by the LHP Coll. [5J, for both the "isovec- 
tor" Dirac ff (obtained through (protonlu^^u — 
d'y tJ 'd\proton)) and the "isovector" radius (with 
different extrapolations), are shown in Fig. [Tj 

3. Phenomenology I: Dispersion Relations 

Dispersion relation (DR) approaches have a 
very long story. Within such a framework, one 
can relate real and imaginary parts of the em 
ff's, once the unitarity (probability conservation) 
and analyticity (causality) properties are imple- 
mented. The analysis based on DR allows one 
to investigate both SL and TL regions, and even 
the unphysical region in the TL region, where a 
possible bound state pp is present, close to the 
threshold, q 2 = 4Mjy. The main ingredient is the 
DR given by (i 
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or a once subtracted version, where Fj(0) (if fi- 
nite) appears 
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The next ingredient is a suitable Ansatz, physi- 
cally motivated, for ImF(t'). For instance, the 
Bonn group [5] has exploited i) the form sug- 
gested by the decay processes p — > irir, </> — ► KK, 
u) — ► pTt, and related continua and ii) the con- 
strains from pQCD. In this case, the fit to the 
experimental data (SL nucleon ff's and TL pro- 
ton data), has 17 free parameters and a total 
X 2 /DOF of 1.8. It should be pointed out that the 
TL neutron calculation represents a prediction of 
the model. In Fig. [2j the TL Nucleon ff's ob- 
tained in the so-called super convergent approach 
(imposed by a proper fall-off of Fi) are shown. 

In Ref. [9] a different approach has been 
followed, by directly applying DR to the ratio 
R(q 2 ) = v P G p E {q 2 )/G p M {q 2 ), both in SL and TL 
regions. A very general form for Im R(q 2 ), with 
six parameters, has been adopted, obtaining a 
X 2 /DOF ~= 1.3. In Fig. El the SL and TL 
ratios are shown (for an analysis of all the ff's 
see Ref. [II])- The approach yields two in- 
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Figure 2. Nucleon ff's for TL momentum trans- 
fers. (After Ref. [8]) 
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Figure 3. R(q 2 ) (see text) vs q 2 . The bands 
represent the error. Green band: retaining in the 
fitting procedure the TL LEAR data [10]. Yellow 
band: the same as the green band, but retaining 
the TL BABAR data [3]. (After Ref. [9]). 



teresting predictions: i) the position of the SL 
zero, that depends upon the TL data chosen 
in the fit, i.e. q 2 z (BABAR) = -10 ± 1 GeV 2 , 
q 2 (LEAR) = -7.9 ±0.7 GeV 2 ; ii) the asymptotic 
limits, extracted from the theoretical expressions 



lim 

q 2 — >±oo 



GW) 



g\M 2 ) 



0.95 ± 0.20 BABAR, 
2.3 ±0.7 LEAR 



4. Phenomenology II: Microscopic Ap- 
proaches 

Relativistic CQ models have been widely 
adopted for obtaining a microscopic description 
of the em properties of the Nucleon. Since '90, it 



was recognized the relevance of dressing the CQ's, 
e.g. by suitable ff's, with a spatial extension of 
the order of 0.4 - 0.5 fm, (a scale different from a 
naive pion-cloud model), in order to obtain a very 
accurate description of the Nucleon and Pion ff's 
in the SL region (see, e.g., [H])- After that, new 
approaches, aimed at a more deep understanding 
of the degrees of freedom beyond the CQ ones, 
have been developed. 

Recently, within the manifestly covariant spec- 
tator model, a very accurate description of the 
Nucleon ff's in the SL region has been achieved 
[13] . The Nucleon, composed of three valence 
CQ's with form factors, has a center of mass mo- 
tion described by a solution of the Dirac equa- 
tion, while the intrinsic dynamics is given by a 
phenomenological S-wave state, consistent with 
the properly symmetrized covariant spectator for- 
malism. The independence of the intrinsic wave 
function upon the direction of the relative mo- 
mentum of an off-mass-shell quark with respect to 
an on-mass-shell diquark allows one to construct 
an internal pure S-wave Nucleon state, and there- 
fore to investigate the Nucleon ff's for a Nucleon 
living in a well defined orbital angular momentum 
state. It is worth noting that the free propaga- 
tion of three quarks is prevented by a suitable 
choice of the intrinsic wave function. The model 
has been further improved by adding a two-pion 
cut for a better description of the ff's at low Q 2 . 
In Fig. [U the Nucleon ff's, obtained through fit- 
ting procedures that take into account the DIS as 
well, are shown. 

Within the Hamiltonian Light- front dynamics, 
where a meaningful Fock expansion can be con- 
sidered, the Pavia group has devised a CQM + 
"meson-cloud" approach (i.e., \qqq) + \qqq',qq)) 
[14] . In particular, the non valence term is de- 
scribed by configurations where a nucleon, or a 
Delta, is coupled to n, p or w, with extended 
coupling constants modeled by phenomenologi- 
cal vertex functions. The momentum dependence 
of the nucleon and Delta wave functions is cho- 
sen according to the Brodsky-Lepage prescription 
(see, e.g. [IS]), while for the mesons a Gaus- 
sian Ansatz is adopted. In a frame where q + = 
<7° + q z =0, only j+ m is necessary, and it is given 
by a sum of one-body currents for bare baryons 




Figure 4. Nucleon ff's, in the SL region, cal- 
culated within the manifestly covariant spectator 
model of Ref. 13J . Dotted line: 4 parameters, ap- 
proximate isospin symmetry in the CQ ff's; short- 
dashed line: 5 parameters, isospin symmetry fully 
broken; long-dashed line: 6 parameters, approxi- 
mate isospin symmetry + two-pion cut; solid line: 
9 parameters, isospin symmetry fully broken + 
two-pion cut. (After Ref. [T3"] ) 



and mesons. Notably, in the chosen frame, j£ m 
is diagonal in the Fock space. After determin- 
ing i) the two parameters in the pion wave func- 
tion by fitting the corresponding ff (for p and 
u> the same momentum dependence is assumed) 
and ii) the cut-off constants in the extended cou- 
pling constants, according to the cloudy bag- 
model, one obtains the remaining 3 parameters 
by fitting the theoretical ff's, to 8 experimental 
values of proton and neutron at low Q 2 (/i p ( n ) , 
Ga(0), G|(0.15 GeV 2 ), G| (M) (0.15 GeV 2 ) and 
G^ (M) (0.45 GeV 2 )). It is worth noting that with 
no 5"-wave, one obtains p p = 2.87, p, n — —1.80, 
r p = 0.877 fm and r 2 n = -0.064 fm 2 . In Ref. 
[14] . a 1% S' component in the Nucleon and Delta 
wave functions has been also considered, in order 
to allow a better description of G^. In Fig. [5] 
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Figure 5. Light-front Nucleon ff's in SL region. 
Long-dashed curve: meson-cloud contribution; 
dotted curve: valence-quark contribution; solid 
curve: the sum of the two contributions; dash- 
dotted curve: the total result with a 1% S'-state 
in (After Ref. [M]). 



the results show that the meson cloud is signifi- 
cant only for Q 2 < 0.5 GeV 2 . 

In Ref. |16j . a Light-front analysis, already ap- 
plied to SL and TL Pion ff [T7], has been ex- 
tended to the Nucleon ff's. The starting point 
is the Mandelstam covariant expression for the 
matrix elements of the em current. It involves: 
i) the Dirac propagators of the quarks, ii) the 
Nucleon Bethe-Salpeter amplitude (BSA), with a 
proper Dirac structure and a scalar function de- 
scribing the four-momenta dependence, iii) the 
quark-photon vertex, to be modeled in order to 
take into account a microscopical Vector Meson 
Dominance (VMD) contribution as well [17. Af- 
ter choosing a frame where q + ^ (in this frame, 
j+ m is not diagonal in the Fock space) , that al- 
lows one to deal with the SL and TL regions on 
the same footing, one can easily integrate on fc~ 
(i.e. projecting onto the Light front), if only the 
Dirac poles are retained. It turns out that the dif- 
ferent terms that naively one should expect from 
physical considerations, can be straightforwardly 



singled out (see [E]). ^ n particular, in the SL re- 
gion, besides the standard triangle contribution 
(or valence term), one has a contribution pro- 
duced by the pair creation (non valence term or 
Z-diagram term). This term involves both a BSA 
with two on-shell quarks, that can be related to 
the valence state, and a BSA with two off-shell 
quarks, leading to a non-valence term in a Fock 
expansion of the Nucleon state. In the TL region, 
there is only the pair contribution. The quark- 
photon vertex contains a bare contribution acting 
in the triangle term (SL only) and a contribu- 
tion acting when the pair production is present, 
namely Zj x a bare term + Zvmd x a VMD 
term, according to the decomposition of the pho- 
ton state in bare, hadronic (and leptonic) contri- 
butions. The momentum dependence of the Nu- 
cleon BSA in the valence sector has been approx- 
imated by a Nucleon wave function a la Brodsky- 
Lepage (i.e., pQCD inspired) [T5] with a quark 
mass m q = 200 MeV and an overall power 7/2. 
Two parameters are present in the wave function, 
one is fixed through the anomalous magnetic mo- 
ments of proton and neutron, (/j^ h = 2.878, Exp. 
2.793, and tf h = -1.859, Exp. -1.913), and the 
other, p, that controls the end-point behavior, by 
the overall fitting procedure (see below). In the 
non valence sector the momentum dependence is 
approximated by 
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where gAB — {m-A m B )/ [P 2 + Mq(A, B)] and r 
is analogous to p. A similar expression holds in 
the TL region. The four adjusted parameters, 
Z{ s = Z[ v = Zy V M , Zy S M , p and r, are obtained 
from the fit of the SL ff's, (even disregarding the 
proton ratio!) with a,x 2 /DOF ~ 1.7. The model, 
illustrated by Figs. [6] and El shows that in the 
SL region the possible zero in G P E is due to the 
cancellation between the triangle term and the 
Z-diagram contribution, generated by the higher 
Fock components of the proton state. Moreover, 
in the TL region one can see that for the pro- 
ton some strength is lacking for q 2 ~ 4.5 and 
- 8 (GeV/c) 2 (as for the Pion Q7]), while the 
neutron remains a challenge. 
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Figure 7. TL Nuclcon ff's. Solid line: full cal- 
culation. Dotted line: bare contribution only. 
Dashed line: full calculation x 2 (a guide for the 
eyes). (After Ref. [15] ) 



Figure 6. SL Nucleon ff's. Solid line: full cal- 
culation. Dotted line: triangle contribution only. 
G D = 1/(1 + |g 2 |/0.71) 2 . (After Ref. [16]) 



5. Conclusions 

Summarizing, Lattice data are becoming more 
and more accurate, but unfortunately only the 
SL region can be studied, and for the time be- 
ing in a limited range of q 2 . Phenomenological 
models, that can address both SL an TL regions, 
appear still very appealing, and could contami- 
nate each other, e.g. by supplying the analytic 
approaches with ingredients microscopically eval- 
uated. In conclusion, one has a set of tools, one 
fundamental and the others phenomenological, 
that can positively interact, allowing to gain a 
more deep insights in the analysis of the Nucleon 
ff's, that represent an intriguing challenge for the 
theory, given the highly non perturbative nature 
of the physical quantities involved in the few-GeV 
region. 
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